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We report the first three-dimensional wavevector analysis of the jellium exchange-correlation (xc) 
surface energy in the random-phase approximation (RPA) . The RPA accurately describes long-range 
xc effects which are challenging for semi-local approximations, since it includes the universal small- 
wavevector behavior derived by Langreth and Perdew. We use these rigorous RPA calculations 
for jellium slabs to test RPA versions of nonempirical semi-local density-functional approximations 
for the xc energy. The local spin density approximation (LSDA) displays cancelling errors in the 
small and intermediate wavevector regions. The PBE GGA improves the analysis for intermedi- 
ate wavevectors, but remains too low for small wavevectors (implying too-low jellium xc surface 
energies). The nonempirical meta-generalized gradient approximation of Tao, Perdew, Staroverov, 
and Scuseria (TPSS meta-GGA) gives a realistic wavevector analysis, even for small wavevectors 
or long-range effects. We also study the effects of slab thickness and of short-range corrections to 
RPA. 

PACS numbers: 71.10.Ca,71.15.Mb,71.45.Gm 



I. INTRODUCTION 

Modern electronic-structure calculations for atoms, 
molecules, and solids usually rely upon Kohn-Sham (KS) 
density- functional theory (DFT),i2 in which only E xc [n], 
the exchange-correlation (xc) energy as a functional of 
electron density, must be approximated. Semi-empirical 
approximations tend to be limited to systems that resem- 
ble those in the fitted data set (typically small molecules), 
but nonempirical ones are constructed to satisfy uni- 
versal constraints and so should have a wider range of 
applicability^ For example, it is expected that a good 
description of chemical reactions at a solid surface re- 
quires a good description of both the molecules and the 
surface. 

Jellium is a simple model of a simple metal, in which 
the valence electrons are neutralized by a uniform posi- 
tive background that extends up to a sharp planar sur- 
face. The apparent success of the simplest density func- 
tional, the local spin density approximation (LSDA), for 
the jellium surface energy^ motivated early interest in 
density functionals and in refinements of the LSDA such 
as the generalized gradient approximation (GGA)^ 

It was therefore a matter of some concern when 
wavefunction-based Fermi HyperNetted-Chain (FHNC) 7 
and fixed-node Diffusion Monte Carlo (DMC)£ calcula- 
tions for jellium slabs (and their extrapolation to infi- 
nite thickness) predicted surface energies considerably 
higher than those obtained in the LSDA. Indeed, DMC 
is usually a gold standard of accuracy. However, it en- 
counters special difficulties for jellium slabs; 9 further- 
more, the large deviations between the available DMC 
and LSDA calculations have been attributed in part to 



inconsistency between the energy of the inhomogeneous 
system and that of the corresponding homogeneous elec- 
tron gas JJiiL Recent approache a 10 ' 11 ^^ 3 -^'^^ have all 
suggested that the actual jellium surface energies are only 
a little higher than those obtained in the LSDA. The jel- 
lium surface-energy story is presented in full detail in 
Ref. 0J. 

In this paper, we perform a detailed analysis of ex- 
change and correlation in jellium slabs, exact at the 
level of the random phase approximation (RPA) , to show 
that the most refined nonempirical density functional, 
the meta-generalized gradient approximation of Tao, 
Perdew, Staroverov, and Scuseria (TPSS meta-GGA) ^1 
can account even for the most long-ranged xc effects at a 
jellium surface. This is a considerable achievement for a 
semi-local functional that is inherently more reliable for 
short-ranged effects than for long-ranged ones. RPA is 
known to be correct at long range; because it has serious 
deficiencies at short-range and, therefore, cannot be com- 
pared to standard versions of the semi-local functionals, 
we use RPA versions of these functionals in this test. 

In order to separate long-range and short-range xc 
effects, we look at the surface contribution to the 
spherically-averaged real-space xc hole, averaged over the 
electron density of the system, and its Fourier transform 
(wavevector analysis). Langreth and Perdew^ showed 
that the exact xc energy of an arbitrary inhomogeneous 
system can be obtained from a three-dimensional (3D) 
Fourier transform of the spherical average of the xc hole 
density, which is a function of a 3D wavevector k. In the 
case of a plane-bounded electron gas, this wavevector- 
dependent spherical average is dominated at long wave- 
lengths (k — ► 0) by the zero-point energy-shift of the 
newly created surface collective oscillations (surface plas- 
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mons) and takes a simple analytical form. This known 
limit has been used to carry out a wavevector interpo- 
lation correction to LSDA, 5 PBE-GGA, 13 and TPSS- 
metaGGA— xc surface energies. The wavevector interpo- 
lation corrections to these functionals were controlle d 13 ! 16 
by using the exact RPA values reported in Ref. HH, and 
led to a consistent set of predicted surface energies.— 

In a DFT context, the RPA is based upon the 
time-dependent Hartree approximation for the density- 
response function but replacing the occupied and unoc- 
cupied single-particle Hartree orbitals and energies by the 
corresponding eigenfunctions and eigenvalues of the KS 
Hamiltonian of DFT. 5 Hence, it describes the exchange 
energy and the long-range part of the correlation energy 
correctly. Essentially exact RPA surface energies were 
evaluated from single-particle LSDA orbitals and ener- 
gies in Ref. [l2|. These calculations provide an accurate 
standard against which approximate density functionals 
(in their RPA versions) can be tested and normed. The 
RPA versions of LSD and GGA were reported in Refs.[l8l 
and[l^, respectively. Because RPA is not self-correlation- 
free, the GGA for RPA correlation is its own meta-GGA. 
The RPA version of the nonempirical TPSS meta-GGA 
was investigated in Ref. Il6l . 

Unless stated otherwise, atomic units are used 
throughout, i.e., e 2 = fi = m e = 1. 



II. THEORETICAL FRAMEWORK 

The exact xc energy, E xc [n], of an arbitrary inhomo- 
geneous system of density n(r) can be obtained from the 
spherical average n xc (r, u) of the coupling-constant aver- 
aged xc hole density n xc (r,r') at r' around an electron 
at r, as follows^ii 6 - 



when 
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E xc [n] = J drn(r)e xc [n](r), 



(1) 



where e[n](r) represents the xc energy per particle at 
point r: 
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£xc[n](r) = 4 / dk du 
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n xc (r,u), (2) 



with 



n xc (r,u) = J dttn xc (r,r'), 



(3) 



d^t being a differential solid angle around the direction 
of u = r' — r. 

The xc surface energy, a xc , is obtained by subtract- 
ing from the xc energy £^ c [n] of a semi- infinite electron 
system the corresponding energy E x ^ ll (n) of a uniform 
electron gas. In a jellium model, in which the electron 
system is translationally invariant in the plane of the sur- 
face, and assuming the surface to be normal to the z-axis, 
one finds 
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7xc(fc)=2— / dzn(z)b xc (k,z), (5) 

7T J-oo 

with kp = (Stt 2 ^) 1 / 3 , n being the background density, 
and 

b xc {k, z) = 4vr / du u 2 [n xc (z, u) - n^ i{ (u)] . 

Jo ku 1 

(6) 

Alternatively, one can introduce Eq. ([3]) into Eq. ([6]) to 
find: 

i r+ k rib r + °° 
MM) = lj I dz'e ik *(*-*') 
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x fi xc (kf,z,z') -n£" lf (fc) 



(7) 



with fcii = \Jk 2 — k 2 , and n xc (fc|i; z, z') and n"" lf (fc) rep- 
resenting Fourier transforms of the coupling-constant av- 
eraged xc hole densities n xc (r, r') and fi™ lf (r, r'), respec- 
tively. At long wavelengths (k — ► 0), one finds the exact 
limit 5 - 



lxc(k) = (u s - ^uj p ) k, 



(8) 



which only depends on the bulk- and surface-plasmon 
energies lj p — (47m) 1 / 2 and lu s — uj p /\/2, and does not 
depend, therefore, on the electron-density profile at the 
surface. 

The spherical average n xc {z,u) entering Eq. ^ can 
be obtained within local or semi-local density-functional 
approximations (such as LSDA, PBE GGA, and TPSS 
meta-GGA) from model o 16 ' 19 ' 21 ' 22 that require knowl- 
edge of the xc hole density n™ lf (it) of a uniform electron 
gas. Alternatively, rigorous calculations of n™ ll (k) and 
the fully nonlocal n xc (k\\; z, z') entering Eq. ([7]) can be 
carried out from knowledge of the A-dependent density- 
response functions Xunif ^-0 an< ^ X (^ll^i z ? z ')j respec- 
tively, defined by adiabatically switching on the e-e inter- 
action via the coupling constant A and by adding, at the 
same time, an external potential so as to maintain the 
true (A = 1) ground-state density in the presence of the 
modified e-e interaction. 23 ! 24 By using the fluctuation- 
dissipation theorem ) 25 ' 26 one finds: 
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and 



rc(k\\' t z,z') = - 



■xn(z) J J 
8(z - z') 



d\ / duj X A (fc|| , ilu; z, z') 



(10) 



With the aim of testing the performance of local and 
semi-local density-functional approximations for the xc 
surface energy, we compare these (local and semi-local) 
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calculations [obtained from Eq. ([6])] to their fully non- 
local counterparts [obtained from Eq. ([7]) with the aid 
of Eqs. ([9]) and (fT0|) ] at the same level of approxima- 
tion, which we choose to be the RPA. On the one 
hand, we evaluate 7a; C (fc) from RPA versions (LSDA- 
RPA, PBE-RPA, and TPSS-RPA) of the local (or semi- 
local) n xc (z,u) entering Eq. ^ based on the RPA xc 
hole density (it) of a uniform electron gas. On the 
other hand, we evaluate 7 xc (fc) from a fully nonlocal ver- 
sion (exact-RPA) of n xc (ku; z, z') entering Eq. ([7]) based 
[by using Eq. (fTU)) ] on the RPA density-response function 
X X (k\\,uj;z,z'). 



III. RESULTS 
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In the calculations presented below, we have consid- 
ered a jellium slab of background thickness a — 2.23 , 
Xf being the Fermi wavelength (\p — 2ir/kp), and back- 
ground density n = [(47r/3)rf] _1 with r s = 2.07. This 
slab corresponds to about four atomic layers of Al(100). 

For the LSD A- RPA calculations, we have obtained 
the RPA xc hole density n™ lf (u) of a uniform elec- 
tron gas either from Eq. ^ or from a non-oscillatory 
parametrization^i For the PBE-RPA and TPSS-RPA 
calculations, we have always used a non-oscillatory 
parametrization of the RPA xc hole density n™ lf (w)j2I 

For the evaluation of the fully nonlocal (exact-RPA) 
7zc(fc) °f Eq. ©, we follow the method described in 
Ref. [l2|- We first assume that n(z) vanishes at a distance 
zq from either jellium edge^ and we expand the single- 
particle wave functions 4>i(z) and the density-response 
function x X (k\\ > UJ \ z ? z ') m sme an d double-cosine Fourier 
representations, respectively. We then perform the in- 
tegrals over the coordinates z and z' analytically, and 
we find an explicit expression for 7 xc (fc) [see Eqs. (IA1|) - 
(|A5|) of the Appendix] in terms of the single-particle en- 
ergies £i and the Fourier coefficients bi s and Xmn{k\\ , lo) of 
the single-particle wave functions <fii(z) and the density- 
response function x A (^|| ; w ; z , z> ), respectively^ We have 
taken all the single-particle wave functions 4>i(z) and en- 
ergies ei to be the LDA eigenfunctions and eigenvalues 
of the KS Hamiltonian of DFT, as obtained by using the 
Perdew-Wang parametrization 1 ^ of the Ceperley-Alder 
xc energy of the homogeneous electron gas.— For the jel- 
lium slab with r s = 2.07 and a = 2.23A^ considered 
here, the exact RPA xc surface energy is found to be 
a xc = 3091 erg/cm 2 , not far from the corresponding RPA 
xc surface energy of a semi-infinite jellium which is known 
to be a xc = 3064 erg/cm 2 J£ 

In Figs. [T]and [2] we have plotted (solid lines) the exact- 
exchange contribution to 7 xc (fc), i.e., 7x(fc), which we 
have obtained from Eqs. (|Aip ~ (|A5p with the quantities 
Xunif(^' w ) an d Xmn(k\\u) replaced by their noninteract- 
ing counterparts xSnif (k> w ) an ^ Xmn(k\\ > w )> respectively. 
Also plotted in these figures are the LSDA, PBE, and 
TPSS calculations of j x (k) that we have obtained by re- 
placing the xc hole densities n xc (z, u) and n"" lf (it) enter- 



FIG. 1: Wavevector analysis ^(k), versus fc/2fcF, of the ex- 
change surface energy of a jellium slab of thickness a = 2.23Af 
and r s = 2.07. Solid and dashed lines represent exact 
and LSDA calculations, respectively. The LSDA calcula- 
tion has been performed either from the actual exchange 
hole density n" mf (u) of a uniform electron gas, which we 
have obtained from Eq. (|9} with Xunif (k, uj) replaced by 
Xunif(ku), (7i SDA ) or from the non-oscillatory parametriza- 
tion of n% nit (u) reported in Ref. |H (< 7 >^ SDA ). The 
area under each curve represents the exchange surface en- 
ergy: cr^ SDA = 2699 erg/cm 2 and a^ ct = 2348 erg/cm 2 . (1 
hartree/bohr 2 = 1.557 x 10 6 erg/cm 2 .) 



ing Eq. © by their corresponding exchange-only coun- 
terparts (dashed lines). 

The LSDA 7^ (k) represented in Fig. Q] has been ob- 
tained by using both the actual exchange hole den- 
sity n" nlf (u) of a uniform electron gas [dashed curve 
labeled 7^ SDA ], which we have obtained from Eq. © 
with x^ nif (/c,w) replaced by Xunif w )j an d the non- 
oscillatory exchange hole density rt™ lf (u) reported in 
Ref. H2 [dashed curve labeled < 7 >^ SDA ]. ^ SBA (k) 
and < 7 >J; SDA (fc) yield, by construction of the non- 
oscillatory exchange hole density n x (u), the same ex- 
change surface energy a x ; they are also almost identical 
in a wide range of low wavevectors, but < 7 >£ SDA (k) is 
considerably less accurate near k = 2kF where the exact 
j x (k) has a kink. This kink is realistic for jellium- like 
systems, but not for atoms and molecules. 

The PBE and TPSS j x (k) represented in Fig. [2] have 
both been obtained by using the non-oscillatory exchange 
hole density n" n (it) reported in Ref. |22j, which yields a 
wrong behavior of j x (k) at large wavevectors. Neverthe- 
less, both the actual exchange hole density n x nll (u) of a 
uniform electron gas (not used in these calculations) and 
the corresponding non-oscillatory exchange hole density 
would yield the same exchange surface energy a x , by con- 
struction, as occurs in the LSDA. 

Figs. Q] and [2] show that while the LSDA j x (k) consid- 
erably overestimates the exact j x (k) at low wavevectors 
(see Fig.[I|, leading to an exchange surface energy <j x that 
is too large, the PBE and TPSS j x (k) are close to the ex- 
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FIG. 2: Wavevector analysis 7^(fc), versus k/2kF, of the ex- 
change surface energy of a jellium slab of thickness a = 2.23Af 
and r s = 2.07. Solid and dashed lines represent exact 
and semi-local (PBE and TPSS) calculations, respectively. 
The semi-local PBE and TPSS calculations have been per- 
formed from the non-oscillatory parametrization of n" nlf (u) 
reported in Ref. I22I The area under each curve repre- 



sents the exchange surface energy: a x 



2155 erg/cm 



2247 erg/cm 2 , and cr° x 



2348 erg/cm 2 



FIG. 3: Wavevector analysis y X c(k), versus k/2kF, of the RPA 
xc surface energy of a jellium slab of thickness a = 2.23Af 
and r s = 2.07. Solid and dashed lines represent exact-RPA 
and LSDA-RPA calculations, respectively. The LSDA cal- 
culation has been performed either from the actual RPA 
xc hole density of Eq. © (7L SDA_RPA ) or from a non- 
oscillatory parametrization of n™ f (u) (< 7 >if DA )-— The 
area under each curve represents the RPA xc surface energy: 
ct lsda-rpa = 3034 erg/cm 2 and a e x T ct ~ RFA = 3091 erg/cm 2 . 



act j x (k) (see Fig. [2]). We note that the peaks of 7 PBE (fc) 
and j x PSS (k) are close to the exact one, a fact which was 
used in the construction of the TPSS exchange hole^ 
and that at larger wavevectors 7 PBE (fc) and 7j PSS (fc) 
nearly coincide, as expected; at lower wavevectors, how- 
ever, the TPSS meta-GGA differs from the PBE GGA, 
leading to a wavevector-dependent "f x (k) that is closer to 
the exact behavior. 

We have also carried out calculations of the exact 7 X (k) 
for increasing values of the background thickness a, and 
we have found that (i) J x (k) is only sensitive to the size 
of the system at wavevectors below the minimum that 
is present in the solid lines of Figs. [1] and[2J and (ii) as 
k — ► the wavevector-dependent j x (k) approaches in the 
semi-infinite limit the profile-independent negative value 
Hx = —1.50 x 10 4 /r 3 i erg/cm 2 ) reported in Refs. d and 

Figures [3J and [4] exhibit the results that we have ob- 
tained for the RPA j xc (k) from Eqs. (|A"Ij) - (|A"5|l (solid 
lines) and within the LSDA-RPA, PBE- RPA, and TPSS- 
RPA (dashed lines). As in the case of the exchange-only 
contributions represented in Figs. [T] and [2l the LSDA 
lxc{k) represented in Fig. [3] has been obtained by using 
both the actual RPA xc hole density fi"" lf (u) [dashed 
line labeled J x i DA ~ RPA ], which we have obtained from 
Eq. ([9]), and a non-oscillatory xc hole density n£" lf (it) 
[dashed line labeled < 7 >^ DA ~ RPA ]; the PBE and 
TPSS 7xc(fc) represented in Fig. 2] have both been ob- 
tained by using a non-oscillatory xc hole density n™ lf (u). 

Figure [3] shows that at short wavelengths with 
k > 2k F the quantities l X c DA ~ RPA {k) (dashed line) 
and ■fxc act ~~ RPA (k) (solid line) nearly coincide, as 



expected^ ' 32 ! 33 The LSDA, however, considerably under- 
estimates Jxc{k) at low wavevectors. This is partially 
compensated by an LSDA "f xc (k) that at intermediate 
wavevectors (around the peak of ^xcik)) is too large. Fig- 
ure |4] shows that the PBE GGA improves j xc (k) at in- 
termediate wavevectors more than at low wavevectors, 
thereby yielding a xc surface energy that is even smaller 
than in the LSDA. From a different perspective* 3 - 4 - the 
too-small PBE surface energy arises from a too-large gra- 
dient coefficient for exchange, but this is repaired by the 
TPSS meta-GGA which uses the proper gradient coeffi- 
cient. Indeed, Fig. [4] clearly shows that the TPSS meta- 
GGA brings improvements over the corresponding PBE 
GGA at both intermediate and small wavevectors, thus 
leading to a wavevector-dependent r y x PSS ~ RPA (k) that 
is very close to r y^. act ~ RPA (k) (solid line) and to an xc 
surface energy a xc that is only slightly lower than its ex- 
act RPA counterpart^ We have obtained similar results 
(not displayed here) for r s = 3, and we have found that 
the errors introduced by the use of noncmpirical semi- 
local density-functional approximations slightly increase 
with r s as expected from the analysis of Ref. [l6|- 

Also represented in Fig. 2] (by a dotted line) is the uni- 
versal (density-profile independent) low-wavevector limit 
of Eq. ([8]). The TPSS-RPA j xc (k) has the virtue that 
not only is it very close to its exact-RPA counterpart 
in the whole range of low and intermediate wavevectors, 
but it imitates the exact low-wavevector limit of Eq. © 
as well. That this limit is also reproduced by the exact- 
RPA 7a; C (fc) of a semi- infinite electron system is shown in 
Fig.O where we have plotted calculations of this quantity 
for increasing values of the background thickness a, from 
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FIG. 4: Wavevector analysis 'yxc(k), fc/2fcF, of the RPA 
xc surface energy of a jellium slab of thickness a = 
2.23Af and r s = 2.07. Solid and dashed lines repre- 
sent exact-RPA and semilocal-RPA (PBE-RPA and TPSS- 
RPA) calculations, respectively. The semilocal PBE-RPA 
and TPSS-RPA calculations have been performed from a non- 
oscillatory parametrization of n),!"' f (u).— The area under each 
curve represents the RPA xc surface energy: cr P ,? B_RPA = 
2959erg/cm 2 , aJ c PSS ~ RPA = 3052erg/cm 2 , and a x T ct ~ RPA = 
3091 erg/cm 2 . The straight dotted line represents the univer- 
sal low- wavevector limit of Eq. (Jsj) . 
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FIG. 5: Wavevector analysis j xo {k), versus k/2kp, of the ex- 
act RPA xc surface energy of jellium slabs of r s = 2.07 and 
various values of the background thickness: a = 8.23Af ('l')i 
a = 2.23Af ('2'), and a = 0.56Af ('3'). The straight solid line 
represents the universal low- wavevector limit of Eq. ((Sj , which 
corresponds to a plane-bonded semi- infinite system (a — > oo). 



a = 0.56Xf to a = 8.23A_f. Furthermore, Fig. [6] shows 
that "f X c(k) is only sensitive to the background thickness 
at very low wavevectors. 

Finally, in order to investigate the impact of short- 
range corrections to the RPA j xc (k), we have plot- 
ted in Fig. [7] the correlation contribution to ^xdk), i.e, 
7c (fc), as obtained in the RPA (LSDA-RPA, TPSS-RPA, 



FIG. 6: Wavevector analysis *y xc (k), versus k/2kp, of the 
exact RPA xc surface energy of jellium slabs of r 3 = 2.07 
and two values of the background thickness: a — 2.23Af 
(solid line) and a — 0.56Af (dashed line). The area un- 
der each curve represents the exact RPA xc surface energy 
^exact-RPA. 3 o9i e rg/cm 2 and 3043 erg/cm 2 , for a = 2.23A F 
and a — 0.56Af, respectively. 
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FIG. 7: Wavevector analysis 7c(fc), versus k/2kF, of the corre- 
lation surface energy of a jellium slab of thickness a = 2.23Af 
and r s = 2.07. Dotted, long-dashed, and dashed-dotted 
lines represent LSDA-RPA, TPSS-RPA, and exact-RPA cal- 
culations, respectively. Short-dashed and solid lines repre- 
sent standard versions of the LSDA and the semi-local TPSS, 
as obtained from an accurate (beyond RPA) non-oscillatory 
parametrization of the correlation hole density n" mf (w) of a 
uniform electron gas.— The area under each curve represents 
the correlation surface energy: ^lsda-rpa _ ggg erg/cm 2 , 



TPSS-RPA 



= 804 erg/cm 2 , a 



exact-RPA 



290 erg/cm , and a, 



TPSS 



756 erg/cm 



= 743erg/cm 



and exact-RPA) and also in standard versions of local 
and semi-local density-functionals (LSDA and TPSS) 
that use an accurate (beyond RPA) non-oscillatory 
parametrization of the correlation hole density n" (w) 
of a uniform electron gas. ?l We observe that in the long- 
wavelength limit (k — * 0), where both LSDA-RPA and 
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standard LSDA exhibit serious deficiencies, both TPSS- 
RPA and the more accurate standard TPSS coincide with 
the exact-RPA. At shorter wavelengths, the standard 
TPSS predicts a substantial correction to its TPSS-RPA 
and exact-RPA counterparts, which is first positive and 
then negative and leads, therefore, to a persistent cancel- 
lation of short-range correlation effects beyond the RPA 
similar to the cancellation that was reported in Ref. [l(] 
in the framework of time-dependent density-functional 
theory and a two-dimensional wavevector analysis of the 
correlation surface energy. 



long-range effects. 
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IV. CONCLUSIONS 



APPENDIX A 



We have reported the first 3D wavevector analysis of 
the jellium xc surface energy in the RPA, and we have 
used this fully nonlocal (esentially exact) RPA calcu- 
lation to test RPA versions of nonempirical semi-local 
density-functional approximations for the xc energy. We 
have tested the first three-rungs of the Jacob's lad- 
der classification of nonempirical density functionals: 36 
LSDA, PBE GGA, and TPSS meta-GGA. 

We have found that while the LSDA displays can- 
celling errors in the small and intermediate wavevector 
regions and the PBE GGA improves the analysis for 
intermediates wavevectors while remaining too low for 
small wavevectors (implying two-low xc surface energies), 
the TPSS meta-GGA yields a realistic wavevector anal- 
ysis even for small wavevectors or long-range effects. We 
have also demonstrated numerically the correctness of 
the LSDA at large wavevectors^^^ (where LSD-RPA, 
TPSS-RPA, and the exact-RPA coincide, as shown in 
Fig. [7]) and the universal low- wavevector behavior derived 
by Langreth and Perdew^ which is nicely reproduced by 
the TPSS meta-GGA. 

We have carried out fully nonlocal RPA calculations 
for increasing values of the background thickness, and 
we have found that the 3D wavevector analysis of the 
xc surface energy is remarkably insensitive to the slab 
thickness except at very long wavelengths (k — ► 0) where 
decreasing the slab thickness reduces the universal slope 
that is dictated by the presence of bulk and surface col- 
lective oscillations. 

Finally, we have found that the TPSS wavevector anal- 
ysis of the correlation surface energy, as obtained from an 
accurate (beyond RPA) non-oscillatory parametrization 
of the xc hole density of a uniform electron gas, provides 
both the exact short- fc limit, where LDA fails badly, and 
the exact large-A; limit, where RPA is wrong. Hence, 
our calculations support the conclusion that the TPSS 
meta-GGA xc density functional accurately describes the 
jellium surface, including not only short-range but also 



Here we give an explicit expression for the wavevector- 
dependent contribution j xc {k) to the xc surface energy 
a xc of a jellium slab of background density n and thick- 
ness a, in terms of the single-particle energies Si and 
the Fourier coefficients bi s and Xmn °f the single-particle 
wave functions <j)i(z) and the density-response function 
X\{k\\;z,z'), respectively^ From Eqs. ©, ©, ©, and 
(|10p and performing the integrals over the coordinates z 
and z' analytically, we find: 



Jxc(k) = — 

TT 



f dk 

/ —r- ot mn (k z )p mn (k\\) 

Jo m=0n=0 

nan*f(k)], (Al) 



where 



2 1 - (-l) m cos(fc z d) 
[Jfc2 - {mir/d) 2 ][k 2 z - (nn/dy 



(A2) 



anct^ 



Pmn(k\\) = -- / dX du;x mn (k\\,iuj) 



~^d? 



X> F - e; )£ G "' G "'< (A3) 



with d = a + 2zq, 



and 



1, for m = 0, 

2, form>l, 



(A4) 



OO OO 



C*W — 9 ^ ^ bish's'i&m.s-s* + 5m, s' 



— s "m,s-\-s' J 



8=1 S' = l 



(A5) 
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